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Abstract 



We here develop a theory on optical absorption in a dilute Bose gas at 
low temperatures. This theory is motivated by the Bogoliubov theory of el- 
ementary excitations for this system, and takes into account explicitly the 
modification of the nature and dispersion of elementary excitations due to 
Bose-Einstein condensation. Our results show important differences from ex- 
isting theories. 

PACS numbers: 03.75.Fi, 32.70.Jz, 32.80.-t 

A remarkable property of a degenerate Bose gas is the modification of the nature of exci- 
tations in the system. This has been elucidated by the work of Feymann and Bogoliubov [|TJ . 
For a classical or non-degenerate Bose gas, elementary excitations are basically quasiparti- 
cles. For a degenerate Bose gas however, the long wavelength elementary excitations are the 
same as density waves of the system, propagating at the sound velocity. They correspond to 
neither addition nor removal of a quasiparticle, but a linear combination of both. At short 
wavelengths these excitations resemble more closely quasiparticles, with the dispersion the 
same as free particles apart from a shift. 

Recently there is much interest in optical excitations in dilute degenerate Bose gas. The 
frequency of absorption has been used as an indication of the formation of Bose-Einstein 
condensation in the system . An optical excitation is a non-trivial process, as one of the 
original atoms (referred to as a-atoms below), all identical before the excitation occurs, is 
converted to another one which is distinguishable from all the others. Moreover a-atoms 
which are not excited internally interact differently with this 'foreign' atom (referred to as 
c-atom below) than among themselves, and must respond by rearrangement of their relative 
motion. The frequencies at which optical absorption occur are therefore different from ujq, 
the value for a single isolated atom. 

Optical absorption in a Bose gas has been considered by Oktel and Levitov || , and also 
by Pethick and Stoof ||. Johnsen and Kavoulakis || considered the case of excitons in 
semiconductors where the initial and final bosons have different effective masses, using an 
approximation equivalent to that used by Oktel and Levitov ||. Let us briefly summarize 
the main results of Ref. most relevant to the discussions below. These authors concen- 
trate on intermediate to high temperatures and ignore the modifications of the nature of 
quasiparticles below the transition temperature. They make the dramatic prediction that, 
for T < T c and in the limit of no momentum transfer, there is actually not one but two 
absorption lines. These lines are located at the frequencies 



u = uj + g ac n - g aa (n + 2n T ) + y 



(1) 



where y obeys the quadratic equation y(y — g ac nr + g aa no) = gl c n o n T- Here g aa and g ac are 
the interaction parameters among the a-atoms and between c- and a-atoms respectively, n 
the number density of the condensate and tlt that of the thermal atoms, n = Uq + the 
total number density. As T — > 0, y — > and — g a an though the weight of the second line 
also approaches zero. For T — > T c , y — > and y — > g ac -n, with the weight of the first line 
vanishing. Note then the second solution of y then yields the line to = u + 2(g ac — g aa )?i. 

Here we would like to re-examine and extend this investigation to low temperatures. In 
particular, we are interested in the effects brought about by the macroscopic occupation 
of the lowest energy state. Note that not only the dispersion is changed for T < T c , but 
also that now the annihilation of an a-atom can correspond both to an absorption and a 
creation of a Bogoliubov excitation. We shall see that properly accounting for these features 
substantially modifies the absorption spectrum. In general it no longer simply consists of 
two lines. In some regimes where the spectrum is approximately two relatively sharp lines, 
the frequencies of these lines differ from those of eq (|T]). 

We then consider a weakly interacting Bose gas. The Hamiltonian of this system is given 

by 

H a = ip a \, a p + 2y 51 a pi a p2 a p3 a p4 (2) 

Here £ p = p 2 /2m is the kinetic energy of a free a-boson of momentum p (we shall drop all 
vector labels since no confusion will arise), and the second sum is over p's with the constraint 
P1+P2 = P3 + P4:- V is the volume. The Hamiltonian of the system containing also a c-atom 
is given by H = H a + H c + H ac . H c describes the free motion of the c-atom 

H c = J2(£p + u o) c I c p ( 3 ) 
p 

where we have neglected the effective mass change in the optical transition, ujo accounts for 
the internal energy difference between the a and c atoms. H ac describes the interaction of 
the c-atom with the rest of the a-atoms. We write it as 

H ac = y- cl,_ q c k ,al,a p ,_ q (4) 

k',q,p' 

The absorption with momentum transfer k (from the external optical perturbation) at 
frequency u is proportional to —1/tt times the imaginary part of the propagator D(u>) whose 
expression in real time is given by 

D(t) = -i < [X(t),X f (0)] > 0(t) (5) 

where X = CfcaJ, + Y^p^o c k+p a j,- 0{t) = 1 if £ > and vanishes otherwise. For here and 
below the angular bracket denotes the equilibrium expectation value at temperature T for 
a system with no c- atoms under the Hamiltonian H. 

As in we are led to consider the equation of motion of the operator X(t) and hence 
CfcaJ and Ck+ P a p . To motivate the procedure below let us consider the operator Ck+ p a p 
and the contribution to its equation of motion from its commutator with H a describing the 
interaction among the a-atoms. The operator Ck+ P behaves simply as a scalar, and the terms 



2 



that are produced, apart from this factor of Cfc +P , are identical with that for the operator 
ap in a BEC of the a- atoms alone. Thus we have a linear combination of terms involving, 

among others, (with Ck+ P implicit) ala a p , alala_ p , and ap,a p iap (p f ^ p). Explicitly, we 
have 

■ 9 . +. / g aa f f n 9aa t t I (a\ 

«^(cfc06) = -Cfc I -^a»o + 2— 2^ Qp/V a o I + ••• ( 6 ) 



i— (c k+p a\) = -c fc+ p I £ p aJ + 2-^aJ(4a + + -^ a - P («o a o)J + ••• ( 7 ) 

where the ellipses represent contributions from the other parts of the Hamiltonian. We shall 
treat the mentioned terms by exactly the same approximation as in Popov's generalization 
of Bogoliubov's theory to finite temperatures ||, i.e., a and al are replaced by scalars ( 
y^I%e TllJ,t ) and a pl a p /a p is replaced by < a p ,a p i > ap where the angular bracket represents 
equilibrium expectation values. Note that thus accordingly the chemical potential fi = 
g a a(no + 2J2 P7 to n p) where N p =< a p a p > is the number of particles at momentum p and 
n p the corresponding density. We are thus led to the conclusion that the operator Cfc +P a_ p 
is generated automatically. This operator was not taken into account in ref [Q. Also, if 
the Hamiltonian were to consist of H a alone, these equations, together with those with ap 
replaced by be diagonalized by the Bogoliubov transformation with the excitation 

operators a linear combination of a p and at_ p (with coefficients identical to those of a system 
of a-atoms only). 

The commutators of c^Oq and Ck+ P a p with H c are simple since a's and a^'s then act like 
scalars. These terms simply describe the free motion of the c-atom. Finally, we consider 
in some detail the contributions from H ac and emphasize another crucial difference of our 
present treatment from that of 0. Consider Cka\H ac . The contributions from the q = 
term in eq (|j) simply gives rise to g ac n times the original operator. For q ^ 0, we have 

Sfc' c kc\i_ q Ck'a\{al q aQ + a\a^ q + J2 P 'yto, q a P > a p'-q) where we have split off explicitly terms 
involving the condensates. Since there is at most one c-atom this expression has a finite 
contribution only when k' — q = k. Thus we have terms of the form Cfc +(? ao(aJao + aja^) 
involving explicitly the condensate. This has the obvious interpretation (see more details 
below): the c-atom can be scattered from k + q to k in two ways in the presence of a 
condensate, either by removing/adding one condensate atom while creating/annihilating 
another at +q/—q. Note that each of these in turn can correspond to emission/annihilation 
of an excitation at +q/—q. Again replacing ao and by scalars implies that one needs to 
take into account the operator Cfc +P a_ p . 

The above leads us to evaluate the equations of motion of the propagators —iy/No < 
[c k (t),X(0)] > 6(t), -i < [c k+p al(t),X(0)) > 0(t), and -% < [c fc+p a_ p (t), X(0)) > 6{t) 
(p ^ 0) . We shall do so under the generalization of the Bogoliubov-Popov approximation 
as explained earlier. Alternatively, we can introduce from the start the new operators a p 
by a p e lllt = u p a p — f p aL p , the same Bogoliubov transformation as in the corresponding 
pure a-atom system ( For simplicity, we shall choose the gauge where u p and v p are all 
real. Note also u p and v p are both even in p.) We then consider the three propagators 
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Tjf = -iVW < [c k (t),X{0)} > 8{t), 7$ = i < [c k+p al(t),X(0)\ > 6{t), and 77$ = 
— i < [cfc +p a_ p (t), X(0)] > 6{t) . This has the advantage that the part due to H a is already 
diagonalized. These propagators are related to the probability amplitude that, at time t, the 
system is in the state with an exciton consisting of, respectively, annihilation of a condensate 
atom at p = with the c-atom at k, annihilation of an excitation of momentum p with the 
c-atom at k + p, and creation of an excitation of momentum — p with the c-atom also at 
k + p. We shall call these 'excitons' type 0, 1 and 2 respectively. The resulting equations of 
motion read, with u' = uo — (uj + g ac n — p). 

(^-^-g^iv^-v^+No (8) 

<y- {£p+k - E p ))rfpl = g a Jy{u p - v p )rif ] + g a Jy{u p r]t ] + v p 7] { k u) ) + u p N p (9) 
(V ~ (W + e p))VpI = 9ac Np y 1 (u p - Vp)i] { k 0) - g ac Np y 1 (vpvl^ + UpVk ] ) - v p (Np + l) 



(10) 

where N p =< aj,a p >= eE ^T_ x is the bosonic distribution function (the number of excita- 



tions, not particles), E p = \JQ + 2£p(gaa,no) the excitation energies under the Bogoliubov- 
Popov approximation, and we have introduced the short hands rjjf 1 = J2 p ^o( u pVpi: ~ v pVpk) 

and r)k = J2p^o( v pVpi: ~ u p%u)- Note that u' is identical with y of eq (0) if one replaces tit 
by S P n p an d uses the mentioned approximation for the chemical potential. 

The interpretation of these equations is clear by an examination of their form. In the 
absence of H ac , the excitons 0, 1, 2 have energies (measured with respect to ojq + p) 
£ p+ fc — E p , and £, p+ k + E p respectively. These different possibilities of excitons are coupled 
together by the interaction g ac . The factors N p and N p + 1 are the bosonic factors associated 
the absorption or emission of an excitation at ±p. Note that even in the absence of the 
Bogoliubov transformation (putting u p = 1 and v p = by hand), rj p I is still finite: it can 
be generated from a type exciton by emission of an excitation to —p, and is allowed so 
long as iVo is macroscopic (note that the definition of 77^ contains the factor v^/Vo- ) These 
type 2 excitons can in turn be converted back to type 0, as explicitly shown in eq (§). These 
processes were not included in the calculations of ref || and ||. 

Let us first check how our results reduce to those of 0] for T > T c . In this regime there 
is no need to treat T] k in a special manner since iVo is no longer macroscopic. Also u p = 1 
and v p = 0. In this case we need only rff k and equation @ becomes equivalent to that 
used in Ref. . For k = there is only one line and the frequency is given by, from eq (0) , 

= g ac n- Since p = 2g aa n, u = 2(g ac - g aa )n, reproducing the results of @] and [g. 

For T < T c we need to consider the three coupled linear equations (|3|-|T0"|). The terms 
on the right hand side involving explicitly g ac represent 'vertex corrections' if the above 
calculations were formulated in terms of Feymann diagrams. They in general cannot be 
ignored. However, it is of interest to pretend that this can be done and examine the results. 
This is a good approximation at very low temperatures for small gaseous parameters ( see 



4 



Appendix). For simplicity we shall discuss only k = below. In this case rj k has a pole at 
u' = 0. T)Q however has a pole whose location depends on p. It is responsible for transitions 
with u' ranging from for p — > to — g aa riQ for large p. ( this is due to the shift between £ p 
and E p as p — > oo). If we ignore the modification of the spectrum at small p's we have in 
total two lines, one at to 1 = and the other at —g aa nQ, with the weight of the second line 
vanishing if T — > 0. This result is similar to that of ref [[J at the same low temperatures 
limit . 

In general we can obtain from eq (|8|-|i0|) three coupled equations involving rif \ rj^ and 
Tjfc . The formulas are lengthy and we shall not show them here. The desired propagator is 
given by D = rj^ + rj^' . 

We shall show the numerical results below. For technical simplicity we shall fix n and 
the (modified) gaseous parameter (noa 3 ) 1 / 2 while vary the temperature. [ We do this to 
avoid solving self-consistent equations for the chemical potential]. In so doing the total 
number of particles n is not fixed when the reduced temperature t = T/g aa no varies ( being 
a function of only t and (noa 3 ) 1 ^ 2 ). However, we checked that in the (low) temperature 
range investigated below the change in the total particle density amounts to less than a few 
percent. Thus our results can still illustrate the semi-quantitative behavior for fixed n. 

An example for g ac / ' g aa = 20 is as shown in Fig. [1]. At very low temperatures t — 0.1 
the absorption is basically a sharp line at dimensionless frequency uj = to' / {g aa n o) — 
except for small wings from absorption (emission) of quasiparticles at frequencies above 
(below) the main absorption line. At increasing temperatures the lower frequency wing 
grows, eventually evolves into something which resembles a sharp line at a frequency which 
decreases with temperature. At the same time the upper line increases in frequency. We 
checked that at higher temperatures than the ones shown, the weight of the lower (upper) 
line decreases (increases). While our calculation cannot be simply extended into higher 
temperatures (see below), we believe that the upper line will eventually evolve to the line 
w = 2(g ac - g aa )n for T > T c . 

The results for a large and negative g ac is as shown in Fig |2]. At very low temperatures 
the results are qualitatively similar to positive g ac . However, at temperatures as low as 
t = 0.4 the results are qualitatively different from those of g ac > 0. The upper line is much 
broader and extends to much higher frequencies. The dependence on the sign of g ac becomes 
even more significant at higher temperatures, as shown in Fig [3] for t = 1.0. Note that for 
the Is — > 2s transition in atomic hydrogen ||, g ac /gaa is likely to be large and negative. The 
implication of this on experiments in traps where the density is non-uniform needs however 
further investigations. 

Though we found basically two lines at intermediate temperatures, we note here that 
the separation between the lines are quantitatively different from those predicted in 0. 
Since the distinctions between quasiparticles and real particles are ignored in H] there is 
an ambiguity in generalizing eq. ([I]) quantitatively, according to whether one uses the total 
number of excitations, or the total number of p ^ particles n'(T), or n'(T) — n'(0) for 
tit- For the first choice we find that, e.g., for g ac / g aa = 20 at t = 1.0, Cj should be given by 
2.49 and —3.1. The third choice yields 3.39 and —3.76, while the second choice produces an 
even bigger separation between the two lines. In all cases the separations between the lines 
are overestimated (c.f. Fig|3|). The disagreement for all choices further worsens at higher 
temperatures. We believe that these differences result from the neglect in p5] of the third 
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channel ( 'type 2 excitons'). 

We have also investigated the dependence of the absorption spectrum on the (modified) 
gaseous parameter (noa 3 ) 1//2 . A smaller (noa 3 ) 1//2 reduces (enhances) the effect due to the 
thermal excitations and the result resembles those of a lower (higher) temperature. 

The treatment here is restricted to dilute gases at relatively low temperatures. At higher 
temperatures one need to worry about more complicated objects involving annihilation and 
creation of multiple Bogoliubov excitations, and the interaction among these various chan- 
nels. The theory for optical absorption will become much more complicated just as a theory 
of the elementary excitations will be at these temperatures. 

This research was supported by the National Science Council of Taiwan under grant 
number 89-2112-M-001-105. This project was motivated during a stay of the author at the 
Aspen Center for Theoretical Physics. The author would like to thank the Center for its 
support. 

Appendix - When solving for 77^ using eq (|TUD, there arises a divergence ~ 
^rr En , , t 1 . jp \ from large p's. This divergence has the same origin as when one studies 
the scattering of two particles interacting via a 5-function potential, and can be cured by 
eliminating g ac in favor of the scattering length. After this divergence is cured, it can be 
shown that as T — ► the vertex corrections are of order (noa 3 ) 1 ^ 2 small. 
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FIG. 1. Absorption spectrum for 
9acl9aa = 20, {n^) 1 / 2 = 0.01, for reduced 
temperatures t shown in the legend. 
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FIG. 2. Absorption spec- 

trum for g ac /g aa = -20, (rtoa 3 ) 1 / 2 = 0.01, 
for reduced temperatures t shown in the leg- 
end. 
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FIG. 3. Absorption spectrum for 
9ac/gaa = ±20, t = 1.0, for (noa 3 ) 1 ^ = .01. 
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